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ABSTRACT 

Given the recently deduced relationship between X-ray temperatures and 
stellar velocity dispersions (the "T-cr relation" ) in an optically complete sample 
of elliptical galaxies ( pavis fc White 1996| ), we demonstrate that L > 
ellipticals contain substantial amounts of dark matter in general. We present 
constraints on the dark matter scale length and on the dark-to-luminous 
mass ratio within the optical half-light radius and within the entire galaxy. 
For example, we find that minimum values of dark matter core radii scale 
as r^m > 4(Ly/3L*) 3//4 /i8 L kpc and that the minimum dark matter mass 
fraction is ^20% within one optical effective radius r e and is ^ 39-85% within 
6r e , depending on the stellar density profile and observed value of P spec - We 
also confirm the prediction of Davis & White (1996) that the dark matter is 
characterized by velocity dispersions that are greater than those of the luminous 
stars: <r^ m m 1.4— 2erJ. The T— a relation implies a nearly constant mass-to-light 
ratio within six half-light radii: M/Ly « 25h 8Q M & /Ly e . This conflicts with the 
simplest extension of CDM theories of large scale structure formation to galactic 
scales; we consider a couple of modifications which can better account for the 
observed T-cx relation. 

Subject headings: galaxies: elliptical — dark matter - X-ray: galaxies 
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1. Background 

There is a strong consensus that spiral galaxies are embedded in massive, non-luminous 
halos, the structure and scaling properties of which have been intensively investigated 



(e.g., [Persic, Salucci, fc Stel 1996|) . Studies of galaxy kinematics and of the density and 
temperature distributions of hot intergalactic gas in galaxy groups ( [Mulchaey et al. 19961 ) 
and clusters (e.g., |Mushotzky 19981) have demonstrated the dominance of dark matter on 



larger scales. In the context of existing theories of the formation of galaxies and large scale 
structure, it is therefore natural to expect non-luminous material within elliptical galaxies 
as well. However, until very recently, the evidence for dark matter in ellipticals has been 
less forthcoming and more controversial. 

Dark matter in ellipticals has been elusive largely due to the lack of internal dynamical 
indicators at optical and radio wavelengths that are as ubiquitous, spatially extensive, 
and unambiguous as HI rotation curves in spirals. Although mergers have endowed some 
early-type galaxies (many classified as SO) with neutral (e.g., [Franx, van Gorkom, fc de 



Zeeuw 1994Q or ionized ( |Bertola et al. 1993|) gas disks, these are rare and may not be 
representative of the population of ellipticals as a whole. Studies of stellar kinematics have, 
until recently, been inadequate in quality and extent to significantly constrain the large 
scale mass distribution in ellipticals in light of the complicating effects of projection and 
anisotropic velocity distributions. 

Most, if not all, giant elliptical galaxies contain extended distributions of hot X-ray 
emitting gas that provide a means for probing the large scale dark matter distribution 

- a means potentially as powerful as rotation curves in spiral galaxies. Since these 
hot gaseous halos should be close to hydrostatic equilibrium (Loewenstein fc Mathews 



1987| ), measurements of the density and temperature distributions from X-ray imaging and 
spectroscopy yield constraints on the mass distribution on the scale and with the accuracy 
and spatial resolution of X-ray temperature profiles. Observations of elliptical galaxies 
by the Einstein Observatory IPC were presented as strong evidence for the presence of 
extensive dark halos in early-type galaxies (|Forman, Jones, fc Tucker 1985[ [Fabian et al.| 



1986| ; |Loewenstein fc Mathews 1987|) ; however, the temperature profiles were of sufficiently 



poor quality to foster a continuing skepticism ( |Fabbiano 1989| , |de Zeeuw fc Franx 1991| ). 



Progress on a number of fronts has been rapid during recent years. In the optical, 
improvements in detector capabilities, data analysis, and modeling have now provided 
robust stellar dynamical constraints on the mass distribution out beyond the optical 
half-light radius r e for several ellipticals QSaglia et al. 199~3| , |Carollo et al. 1995| , |Rix et| 



al. 1997 , Gerhard et al. 1998j ), while gravitational lensing observations are providing the 



first evidence for dark matter in intermediate redshift galaxies flKochanek 1995| ; jBrainerd 
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Blandford, & Smail 1996; [Griffiths et al. 1996 ). X-ray observations using the ROSAT 
and ASCA satellites have improved the accuracy, spatial resolution, and extent of derived 
hot gas density and temperature profiles. Of particular interest is the recent study of 
ASCA observations of ~ 30 early-type galaxies by Matsushita (1997) that shows that 
gas temperature profiles are relatively flat out as far as can be reliably measured (5-20', 
depending on the intrinsic luminosity and galaxy distance). The existence of isothermal 
0.5-1 keV gaseous halos extending, in some cases, to well beyond 10r e is strong evidence for 
the presence of extended dark matter in elliptical galaxies. 

The case for dark matter halos in a number of elliptical galaxies is now overwhelming, 
and the following pattern seems to be emerging from the latest results on individual 
galaxies: dark matter comprises a significant but not dominant fraction of the total mass 



within r e flSaglia et al. 1993; 


De Paolis, Ingrosso, & Strafella 1995 


et al. 1998|; but, see Mushotzky et al. 1994 


for an exceptional case 



becomes increasingly dominant at large radii: M/r ;«5x 1O 12 M kpc -1 to within a factor 
of two for r > 30 kpc (Mushotzky et al. 1994] , [Irwin fc Sarazin 1996| , [Buote fc Canizarel] 
|1997| ), corresponding to mass-to- light ratios > 100M o /L Bq measured at (or extrapolated 
to) ~ 100 kpc ( [Bahcall, Lubin, fc Dorman 1995| ; [Griffiths et al. 19960 . 



2. Motivation and Overview 

While there is now little doubt that some elliptical galaxies do contain dark matter, a 
number of broader issues remain to be addressed. Among the most important with respect 
to our understanding of galaxy formation is whether dark matter in ellipticals is ubiquitous, 
and how the relative scaling of dark and luminous matter varies — both within galaxies as 
a function of radius, and between galaxies with different optical properties. In this paper we 
use the recently derived collection of X-ray temperatures in a complete, optically selected 
sample from Davis & White (1996; hereafter DW) to constrain the average properties of 
dark halos in the population of L > L* elliptical galaxies. 

DW analyzed X-ray spectra from 42 of the 43 optically brightest elliptical galaxies, 
using ROSAT PSPC data if available and Einstein IPC data otherwise [i.e. for four 
galaxies). See DW and White & Davis (1998) for details about the sample and data 
analysis. For the purposes of this work, we use only the mean relationship between X-ray 
temperatures, (T), and central projected optical velocity dispersions, (a) (hereafter, the 
"T-cr relation") — (T) oc (cr) 1 ' 45 — or equivalent X-ray/optical correlations. More precisely, 
we consider the variation of the observable f3 spe c with optical luminosity Ly (Figure 1), 
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where 

Hspec — k(T) ' 

and /im p is the mean mass per particle. If dark matter dominates the gravitational 
potential on large scales, then (T) is essentially a measure of the dark matter content 
within the extraction radii used by DW (6r e , where r e is the optical effective radius). 
Since the observed scaling relations (i.e. the fundamental plane) for elliptical galaxies 
provide a link between (a) and the global luminosity, /3 spec is an excellent diagnostic of the 
dark-to-luminous matter ratio within the optical radius. 

Temperature profiles can provide additional dark matter constraints, but are 
measurable for only the X-ray brightest galaxies in the sample. Likewise, extended stellar 
velocity dispersion profiles are useful, but data of sufficient detail to be unambiguously 
interpreted are not generally available. On the other hand, central velocity dispersions have 
been accurately measured for the entire sample; central velocity dispersions are known to 
correlate with other optical properties and are less likely to suffer from the uncertainties 
caused by anisotropic stellar orbits. 

Our goal in this paper is to make general inferences about the presence of dark matter 
halos and their mean systematic variation with optical luminosity. Toward that end, we 
utilize a set of physically and observationally well-motivated mass models and scaling 
relations in order to reproduce the T-a relation. Although we address the scatter in the 
T-a relation, a detailed examination of individual galaxies is beyond the scope of this paper 
and is deferred to future investigation. 

Our results can be considered in three stages. First, we demonstrate how the observed 
range of /3 spcc necessarily requires the existence of dark halos in elliptical galaxies. Second, 
we derive conservative lower limits on the dark matter scale length and dark-to-luminous 
mass fraction, and discuss how the dark matter content must scale with optical luminosity 
so as to reproduce the observed T-a relation. Finally, we discuss possible implications of 
our results for cosmology and galaxy formation, by embedding our models within particular 
scenarios for large scale structure formation. 



3. Modeling 

As stated above, we use the observable P spec , defined in equation (1), to constrain 
the dark matter distribution in elliptical galaxies. This parameter is the ratio of suitable 
averages (to be defined below) of the square of the stellar velocity dispersion and the 
gas temperature. Therefore, we must derive the expected distributions of stellar velocity 
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dispersions a and gas temperatures T, given specified total gravitational mass distributions. 



3.1. The Stellar Component 

Recently published Hubble Space Telescope (HST) observations of the centers of 
elliptical galaxies show that their surface brightness profiles can generally be characterized 
by the function 



I{R) = 2 (6-c)/a J(rbr) ( M£ 



n n (c—b)/a 



(2) 



where R is the projected radius and r br is a break radius such that / oc R c for R <C r br 
and I oc R~ b for i? r br , while a characterizes the sharpness of the break ([Faber et al. 



|1997|) . In order to accurately calculate the central stellar velocity dispersion, this observed 
departure from distributions with "analytic cores" must be accounted for. In the Appendix 
we demonstrate that, for a = 2 and c < 1, equation (2) provides a good approximation to 
the projection of a space density profile oc r~ d (c < d < c + 1) for r <C r br and oc r~ b ~ l for 
r 3> r br . To derive the luminosity density distribution l(r), we numerically deproject the 
surface brightness profile: 



U r ) = dz(l-x 2 z 2 )- y z(l + z a ) ( )/a (cz a + b), (3) 

7rr br Jo 

where x = r/r br , and the change of variables z = r br /-R has been made. 

The DW sample is magnitude-limited, and we are primarily interested in the most 
luminous ellipticals, since they have the most accurate X-ray temperatures (see §3.3). We 
adopt the typical values for My < —22 galaxies (adopting H = 80 km s -1 Mpc -1 ) of 
a = 2, b = 1.44, and c = 1/10 ( [b'aber et al. 199"7| ). For the nine galaxies in common with 
DW (excluding the highly flattened NGC 4697), the mean surface brightness parameters 
and variances from Faber et al. (1997) are as follows: a = 1.88 ± 0.51, b = 1.38 ± 0.16, 
c = 0.096 ± 0.076. The resulting stellar density profile diverges slightly faster than r _1 
for r C r|, r . The asymptotic surface brightness slope b is flatter than in conventional 
surface brightness parameterizations and the integrated luminosity L(r) does not converge. 
Therefore, we truncate the optically luminous distribution at radius r max , determined so 
that 

L C ore = 7rr^ r /(r br ) = 0.012L(r max ); (4) 

consequently r max = 200r br . 
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As an alternative stellar model, we consider the density profile from Hernquist (1990), 

where L is the total luminosity and rn cin = 0.45r e (as adopted here) provides the best match 
to a de Vaucouleurs (deV) model with effective radius r e ( [Hernquist 1990|) . An advantage 
of this model is that the density and mass profiles are analytic. We shall subsequently refer 
to the stellar distribution of equation (3) as the "HST" model, and that of equation (5) as 
the "Hernquist" model. A comparison of the luminosity density (in units of the average 
density, 3L/47rr max 3 ) as a function of normalized radius (r/r max ) in these models with the 
corresponding deV model is shown in Figure 2, where r max = 200rb r = 40rH ern /3 = 6r e . All 
three models produce similar slopes at very small r; however, deV and Hernquist models 
fail to reproduce the flattening at ~ tv observed by HST. On the other hand, while ~ 80% 
of the total luminosity is contained within m 0.7r max « 4r e in all three models, the HST 
model concentrates the entire remaining « 20% of the mass between w 4r e and 6r e as 
opposed to spreading it out from « 4r e to oo. In order to accurately characterize the 
observed central stellar properties as required by the definition of j3 spec without introducing 
additional parameters (i.e. an additional break radius and a distinct asymptotic form for 
the density profile), such a redistribution of the stellar light is necessary. If the light profile 
does indeed steepen outside of the HST field-of-view, then the only effect on our results 
for the HST model is that, by overestimating the stellar mass at large radii, we would 
underestimate the dark halo mass required to reproduce any observed value of f3 spec . 

The stellar velocity dispersion is derived by solving the Jeans equation in the form 

d ( 2\ 2 GM{<r) . r . 

-j- \P*°r = p*cr r 5 p*, (6) 

ar v ' r r z 

where p* is the stellar density distribution obtained from multiplying the luminosity density 
(see equations 3 and 5, above) by a constant stellar mass-to-light ratio; M(< r) is the 
total mass within radius r (the dark matter component of which is described in the next 
section); A is the velocity dispersion anisotropy parameter (A = l — a t 2 /a r 2 ); and a r (at) 
is the radial (tangential) component of the velocity dispersion. We make the reasonable 
assumption that the velocity dispersion becomes isotropic at small r but radial at large r 
( |Ciotti fc Pclligrini 1992 ), and adopt the functional form 



y*2 



A = =23^5' W 



r z + s 

where the transitional scale radius s is a free parameter of our models. 
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3.2. The Dark Matter Component 

The dark matter mass distribution is assumed to follow the "universal" density profile 
fllNavarro, Frenk, fc White 1997| ), 



Pdm OC ( J I 1 + 

> '"dm / V '"dm 



that, when integrated, yields 

M dm oc f(r/r dm ) = 



, r + r dm 
In- 



^dm r + r dn 



(9) 



It is not clear how appropriate equations (8) and (9) are for elliptical galaxy dark 
halos. On one hand, there is evidence for cores in spiral galaxy dark matter distributions 
(e.g., Persic, Salucci, fc Stel 19961) ; on the other, dissipational collapse of the baryonic 



component in ellipticals can increase the dark matter concentration. Moreover, models 
with cores cannot generally be consistently coupled with stellar mass distributions as cuspy 
as those often observed in elliptical galaxies ( |Ciotti fc Pclligrini 1992| ). We have verified 
the existence of self-consistent two-component mass models with dark matter following 
equation (8) and Hernquist stellar models (equation 7), in the parameter range of interest, 
according to the criterion of Ciotti and Pelligrini (1992). 

From the virial theorem it follows that the integrated mass can be robustly estimated 
from the mean temperature. Constraints on the dark matter concentration can also be 
placed from consideration of the T-a relation; however, any inferences about dark matter 
at small or large radii are contingent on the corresponding assumed asymptotic dark matter 
density slopes. Our choice of equation (8) for the dark matter distribution does allow us 
to connect our results with predictions of dark halo structure from numerical structure 
formation simulations. However, the true shape of the dark matter spatial distribution can 
only be constrained from temperature and stellar velocity dispersion profiles of sufficient 
extent and resolution. 



3.3. The Hot Gas Component 

The hot, X-ray emitting gas is assumed to be in hydrostatic equilibrium so that 

± (p^kT\ = GM(< r) ^_ 

The mass in hot gas is not significant compared to the stellar mass, so its contribution to 
the gravitational potential is neglected. Since we are interested in the average properties of 
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the DW sample, and because spatial analysis is not available for most of the sample, we 
assume that p gas 2 oc P*', this has been found to be a good approximation for X-ray bright 
galaxies, and is expected if heating is balanced by cooling locally ( [Loewenstein fc Mathews 
1987] ) . Both assumptions, hydrostatic equilibrium and p gas 2 cx p*, are more likely to break 



down for galaxies with low X-ray-to-optical luminosity ratios, where a partial or transient 
outflow might obtain; however, the T-o relation for these relatively X-ray faint systems 
does not differ, within the statistics of the sample, from that of the sample as a whole. 



core 



3.4. Scaling Relations and Model Solutions 

A significant unifying simplification for the HST stellar models follows from the 
observations that the ratios of core mass to total stellar mass (M core /M*(r max ), where M< 
is L core [equation 4] multiplied by the stellar mass-to-light ratio), and break radius r br to 
r e (the "true" r e as observed, not the half-light radius in our models) are approximately 
constant: M core ss 0.012M*(r max ) and r Dr ~ ?"e/33 ( [Faber et al. 1997]) . Our adoption of 
these scaling relations, along with the empirical relation M^/Ly oc Ly 1//4 , insures that 
our model galaxies are consistent with the fundamental plane (FP) as long as the core is 
baryon-dominated, since it follows that 

a * M^l o.Ml^Lyo, a c 8 /V e 4 / 5 , (11) 

^br ' e 

where o c is the central velocity dispersion.^ In order to provide us with a one-parameter 
family of models (the quality and quantity of the X-ray data is insufficient to attempt to 
derive a "T-cr-r e " or equivalent FP-like relation), we take a "Faber- Jackson cut" through 
the FP such that Ly oc a 4 — * r e oc Ly 4 , normalized so that r e = 6.3 kpc (and o = 250 km 
s" 1 ) at My = -22 for H = 80 km s^ 1 Mpc" 1 . The root -mean-square deviation of galaxies 
in the DW sample from this projection of the FP is only 34 km s" 1 (13.5% of the mean). 
The Hernquist models (equation 5) likewise adhere to the FP, since we universally adopt 

^Hern = 0.45r e . 

Equations (6) and (10) are solved in dimensionless form. There are two fundamental 
(luminosity-dependent) dimensional quantities that we denote as r* (chosen to be r D r for 



2 We note at this point that we do not distinguish between "dark" ("luminous") and "non-baryonic" 
("baryonic") matter. Elliptical galaxies may very well have an additional baryonic dark component in the 
form of stellar remnants that may account for the variation with luminosity of the central mass-to-light 
ratio ( Zepf fc Silk 1996| ). Any such contribution is subsumed into our definition of the stellar ("luminous") 
component through M^/Ly- 
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the HST model or r He rn for the Hernquist model) and M (M core for the HST model or total 
stellar mass M* for the Hernquist model). Thus, adopting the above scaling relations, the 
length-, mass-, density-, velocity-, and temperature-scaling relations are as follows: 

r* oc X 3/ % -\ (12) 

M oc A 5/ V*/i 8( T\ (13) 

Po = A M \, oc A"V*^80 2 , (14) 



^GM \ 1/2 



^ =\ — ) * A 1/4 /i* 1/2 , (15) 



and 

kT = jj,m p a 2 oc A 1 / 2 //*. (16) 

Here /iso is the Hubble constant in units of 80 km s _1 Mpc -1 , A = Ly / Lq with 
Lo = 5.2 x 10 10 /?,8o~ 2 Lvq, and /x* is defined such that the stellar mass-to-light ratio 
M^jLy = lO/z*M /Ly Q for L = L and h g0 = 1. L h 80 2 corresponds to My = —22, i.e. 

The dark matter mass distribution is scaled accordingly, i.e. 

M dm M,( 

^max I /i r-,\ 

14 =a /(wA) M ' 1 j 

where a is the ratio of dark to stellar mass at the HST stellar model cutoff radius, i.e. 

M dm 



a = 



18) 



and 5 = rd m /^* is the dimensionless dark matter scale length. In the above, x = r/r*, 
x max = r max /r„, and the function / is defined in equation (9). 

To determine the average values of T and a, the hydrostatic and Jeans equations 



are integrated inward from a zero-pressure boundary at r = oo (see |Loewenstein 1994]) . 
Although the stellar mass has a cutoff (r max ) for the HST model, the gas and dark matter 
halos are likely to extend beyond this and are not truncated in the models. This results in 
a pressure contribution to the gas virial temperature within r max that leads to lower values 
of f3 S pec for a given a, and is therefore conservative in the sense of minimizing the required 
amount of dark matter. 

The average quantities that constitute the observable /3 spec are calculated as follows: 
we define (a) as the average projected stellar velocity dispersion within a circular aperture 
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of size r br , and (T) as the emission-weighted projected gas temperature within r max (the 
same aperture over which the temperature is measured in DW, corresponding to the radius 
enclosing m 90% of the optical light). The value of /5 spec is robust to the precise choice of 
these radii since r max r e r^. The aperture-averaged, projected velocity dispersion for 
all central radii > 0.1rb r (> 0.3rt, r ) is within 10% of its value at r^ r for the HST (Hernquist) 
model, assuming isotropic orbits and that the central potential is dominated by the stellar 
component. By averaging over constant multiples of the fundamental stellar scale length 
(r*), we assure that our computed values of /9 spcc depend only on the relative (to stellar) 
dimensionless dark matter and velocity dispersion anisotropy parameters. These parameters 
are the following: the ratio of dark to luminous matter scale length, 5 = rd m /">"*; the mass 
ratio M^m/M* within r max , a; and the dimensionless anisotropy scale length, s/r*. The only 
additional model parameters are the stellar distribution slopes. For a given triplet (s/r*, 
[Md m /M*] rmax , rdm/^*), a general solution for the dimensionless gas temperature and stellar 
velocity dispersion radial distributions follows that can be scaled to physical units for any 
h 80 and Ly using equations (12-16); /i* is determined by normalizing to (a) = 250 km s _1 
for My = —22 and h$ — 1; note that /i* oc (er / (o~)) 2 . 

4. Results 

In this section we consider the dependences of (3 spec on s/r*, rdm/?"*, and [Md m /M*] rmax 
for both HST and Hernquist models. While the latter provides a more accurate global 
representation of elliptical galaxy stellar distributions, the former is more accurate around 
and inside the break radius where (a) is calculated. For models that are dark matter 
dominated beyond several r e , the redistribution of the stellar mass at large radii in the HST 
model is irrelevant. 



4.1. Variation with Anisotropy Scale & the Necessity for Dark Matter 

In the absence of dark matter, the only free parameter in our models is the anisotropy 
length scale s; recall that this is defined such that orbits become isotropic for r <^ s and 
radial for r ^> s. Therefore, if stars dominate the gravitational potential within r max , and if 
the shape of the velocity ellipsoid does not vary systematically with Ly, then p spec should 
be independent of L v and (a) (i.e., (T) oc (<?) 2 ), contrary to the observed T-a relation 
(DW). Moreover, as we now show, the typical observed value of (3 spcc ~ 0.5 cannot be 
reproduced in the absence of dark matter. 
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The solid and dotted curves in Figure 3 show the variation of (3 spec with the 
dimensionless anisotropy scale length s/r^ T for, respectively, the HST and Hernquist stellar 
models without dark matter. Evidently, (3 spec falls with increasing s, but is essentially 
constant for s > 10-20rb r (~ r e /3-2r e /3), since (a) is an average over a much smaller radius. 
The asymptotic values of f3 spec are 1.2 and 0.75 for the the HST and Hernquist stellar 
models, respectively. The lower value for the Hernquist model is the result of a temperature 
maximum at several tv that substantially contributes to the emission-averaged temperature 
because of the high stellar (and, by assumption, gas) density there. As dark matter is 
introduced, the temperature profile becomes more nearly isothermal, and differences in P spec 
between HST and Hernquist models lessen (see next paragraph, and Figure 3). Varying the 
slope parameters of the HST models (equation 2) does not result in any values of (3 spec less 
than 0.75. 

The values of /3 spcc derived without dark matter lie well above what is observed in a 
typical elliptical galaxy (in fact, no galaxy in the DW sample has a best-fit value of P spec 
greater than 1; see Figure 1). One of our major conclusions is already apparent. While a 
typical bright elliptical galaxy has /3 S pec ~ 0.5, even under the most extreme assumptions, 
(3 sp ec is never less than pa 0.75 without the presence of an extended dark halo. While 10 of 
the 30 galaxies in the DW sample have upper limits to j3 spec greater than 0.75, only 3 have 
upper limits greater than 1.2 (see Figure 1). Of the 16 galaxies where (3 spcc is determined to 
better than 20%, only 1 has an upper limit greater than 0.75. 

The dash-dotted and dashed curves show the corresponding models containing dark 
matter halos with dark matter scale length r^m set equal to r e {i.e. r^m = 33rt, r = 2.2rHem), 
and with [M dm = 3M*] rmax chosen so that (3 spcc pa 0.5 for s > r e . The non-luminous 
fraction of the total mass within r e is pa 0.6 and pa 0.5 for the HST and Hernquist models, 
respectively. 

From here on we assume isotropic stellar orbits. As shown in Figure 3, for the practical 
purposes of calculating /3 spec , this is equivalent to assuming that the orbits are mostly 
isotropic inside ~ r e /2. Since f3 spec increases with decreasing s/r br , more dark matter is 
required to obtain a given value of f3 spec if this assumption is discarded — once again we 
adopt the most conservative assumption as regards the required amount of dark matter. 

4.2. Variation with Dark Matter Scale 

Figure 4a shows [3 spcc as a function of the ratio of dark matter scale length to stellar 
break radius, r dm /r br , for the HST (solid curve) and Hernquist (dotted curve) models with 
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[M dm = 3M*] rmax . (For the Hernquist model, the scaling relation r br = r e /33 = r Hc rn/15 has 
been utilized.) As r dm decreases below ~ 10r br , dark matter plays an increasingly important 
role in determining (a), so j3 spcc rises. For fixed [M dm /Mj rmax , /3 spec is nearly independent of 
^dm/^br provided that r dm > r e . It is clear that measurements of /3 spec can only place lower 
limits on r dm . 

There is a fair amount of freedom in how the dark matter can be distributed and 
still produce f3 spcc ~ 0.5. This is illustrated in Figure 4b, which displays the variation in 
baryon fraction at r e (solid curve for the HST stellar model, dotted curve for the Hernquist 
model) and at r br (dashed curve for the HST stellar model, dot-dashed curve for the 
Hernquist model). /3 spec changes very slowly with r dm /r br for r dm /r br > 20, while the baryon 
fraction inside r br varies from ps 0.73 and ps 0.66 for the HST and Hernquist stellar models, 
respectively, to pa 0.98 for both models as r dm becomes very large. Meanwhile, the baryon 
fraction inside r e = 33r br increases from pa 0.31 (pa 0.43) if r dm /r br = 20, to pa 0.75 (pa 0.84) 
for the HST (Hernquist) model as r dm becomes large compared to the stellar cutoff radius 

'"max • 

4.3. Variations with Dark-to-Luminous Mass Fraction 

Figure 5 shows (3 spcc as a function of [M dm /M*] rmax for HST and Hernquist stellar 
models with r dm set equal to 0.5, 1.0, and 2.0r e . As discussed above, (3 spcc — * 1.2 (0.75) 
as M dm — > for the HST (Hernquist) model. The curves pass through j3 spcc pa 0.5 at 
[M dm ^ 3M*] rinax if r dm = r e or 2r e , while [M dm pa 4(6)M*] rmax is required if r dm = r e /2 
for HST (Hernquist) models. As [M dm /M^] rmax becomes sufficiently large for fixed r dm , 
dark matter comes to dominate the mass even at small radii (i.e within r br where (a) is 
calculated) and {3 spec becomes independent of a. That is, in order to maintain consistency 
with the FP and the Faber- Jackson relation, the physical mass must remain constant and 
the increase in dark-to-luminous mass ratio a is compensated for by decreasing the stellar 
mass-to-light ratio (i.e., /i* in equations 13-16) accordingly. 

4.4. Limits on Dark Matter Parameters 

We have shown (Section 4.1) that the lower limits on f3 spec for the HST and Hernquist 
stellar models without dark matter are 0.75 and 1.2, respectively. To explain observed 
values of f3 spec that are less than this requires the addition of dark matter halos with 
sufficient mass to raise the value of (T), but with a large enough scale length to leave (a) 
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(which is averaged over a much smaller scale) relatively unaltered. By calculating solutions 
to equations (6) and (10) for the entire ([M dm /M*] rmax , r dm /r,,,) parameter space, we have 
derived limits on the amount and scale length of any dark halo able to reproduce a given 
observed value of /3 spec - 

As discussed in Section 4.3 and shown in Figure 5, for a fixed value of r dm /rb r there is 
a minimum value of (3 spec that is obtained when a becomes sufficiently large that both (a) 
and (T) are determined by the non-luminous mass component. Since this minimum /3 spec 
increases as r dm /rb r decreases, there is a minimum value of r dm /rb r required to explain any 
observed f3 spec . Figure 6a shows these minimum values of r dm /rb r for the HST (solid curve) 
and Hernquist (dotted curve) models. The minimum value of r dm for a Ly rs L « 3L* 
galaxy with j3 spec ~ 0.5 is ~ 2/i 80 _1 kpc (« r e /3); the dimensional value scales as Ly 3//4 . 
These extreme models are dark matter dominated, even inside the break radius rb r> and 
have very low stellar mass-to-light ratios. Requiring models to have M*/Ly > 5M /Ly for 
Ly = L and h 80 = 1 (or, equivalently, that ^60% of the mass within be stellar) roughly 
doubles the minimum values of r dm /rb r for 0.4 < (3 spec < 0.6. This reasonable requirement 
then implies that r dm > A(L V /SL^^h^ kpc. 

The maximum baryon fractions within r^ r , r e , and r max for any given (3 spec . are shown 
in Figure 6b. The lower limits on the dark matter fractions in models required to explain 
Aspec > 0.4 are less than 5% within r br , and less than 20% within r e — non-luminous halos 
could be quite inconspicuous unless data well beyond r e are considered. The maximum 
baryon fractions at r max for the HST (Hernquist) model are 0.43 (0.15), 0.53 (0.27), and 
0.61 (0.45) to produce p spec = 0.4, 0.5, and 0.6, respectively. These models have very large 
dark matter scale lengths (r dm 3> w) and are analogous to "maximum disk" models of 
spiral galaxies in that they maximize the contribution of stars to the mass-to-light ratio. 



4.5. Scaling of Dark Matter Content with Optical Luminosity 

In order to investigate the scaling characteristics of dark halos with optical luminosity, 
we consider the equivalent of the T-a relation in the j3 spec -Ly plane. Correlations of f3 spec 
and log(/? spec ) with log(Ly) are derived using the DW data set, blue magnitudes and B — V 
colors from the RC3 catalog flde Vaucouleurs, C, et al. 1991Q , and velocity distances 



from Faber et al. (1989). (The (3 spec -\og{Ly) data and best-fit correlation are re-plotted 
in Figure 7b.) Galaxies in the centers of cluster cooling flows and those dominated by 
nuclear point source emission are excluded. The log(/3 spec )-log(Ly) correlation, equivalent 
to P spcc ~ 0.58A ' 13 where A is the visual luminosity in units of L = 5.2 x 10 10 /i 8 o~ 2 Ly Q , is 
essentially identical to the convolution of the T-a correlation and Faber- Jackson relation 
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(see Section 3.4). Nevertheless, requiring that models conform to the observations in 
this form makes use of the tilt of the FP for scaling of the the stellar component, thus 
circumventing concerns about the thickness of the projection of the FP onto the Ly-a 
plane. Note that the galaxies considered span 0.15 < A < 3, although most of the sample 
lies within the range 0.5 < A < 2. 

Given our homologous elliptical galaxy models, which are based on the scaling 
properties of the FP, the observed variation in f3 spec implies that there is a systematic 
variation of relative dark matter properties which breaks the self-similarity for galaxies 
of different luminosity. The sense of the observed T-a relation could also be explained 
by a fine-tuned systematic variation of the velocity dispersion anisotropy (having more 
radial orbits for more luminous galaxies), or if more luminous galaxies have cuspier 
stellar distributions (larger values of c in equation 2). However, the latter is contrary 
to observations; and as we have shown, the observed values of [3 spcc are too low to be 
explained without the introduction of dark matter. It is possible to explain the fl spec -Ly 
trend, assuming the dark-to-luminous mass ratio within r max is constant, if the ratio of 
dark matter to stellar break (or, equivalently, effective) radius, r^ m /r^ T decreases with 
increasing galaxy luminosity Ly. However, since (3 spec is flat for r dm /r^ T > 10 and steeply 
rises for rdm/^br < 10 (Figure 4a) this requires fine-tuning. Moreover, as we discuss shortly, 
7"dm/?"br < 10 is not expected. 

More naturally, the increase in (3 spec with Ly could be explained if [M<j m /M*] rmax 
decreases with increasing Ly (that is, if more luminous galaxies are less dark matter 
dominated) and rdm/^br is relatively large (see Figure 5). We now discuss the plausibility of 
such a correlation within the context of the hierarchical clustering formalism of structure 
formation, simulations of which produced the dark matter distribution expressed by 
equations (8) and (9). This formalism provides specific predictions for r dm as a function of 
virial mass (defined here as the mass within a radius where the mean dark matter density 
is 200 times the average density in the universe) for a given cosmology and power spectrum 
of initial fluctuations. 

In order to extend our models to the scale of the virial radius and mass, we transform 
the dark matter parameters rdm/'"* and [Md m /M*] rmax that are directly constrained by 
observations into the following pair of global parameters: the dark matter concentration 
parameter 

C = r 20 o/r dm , (19) 
and the global dark-to-luminous mass ratio 



^virial = [Mim/Mi<] r200 , 



(20) 
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where r 20 o refers to the radius within which the overdensity of the dark matter (only) is 
200. (M dm (r 20 o) and r 200 refer to the original mass of the dark halo, but this may have been 
subsequently reduced in rich environments through tidal stripping.) Note that for C — 10 
and ctvMai = 10, r 20 o = 320/iscT 1 kpc, r dm /r br = 170, and [M dm /M*] rmax = 1.6 for A = 1 
(L « 3L*). 

In bottom-up hierarchical clustering theories of structure formation, lower mass halos 
are more concentrated — a reflection of the higher density in the universe at the earlier 
epoch of their formation ([Navarro, Frenk, fc White 1997| ). We parameterize this dependence 
with the function 

C = %1- ( 21 ) 

for the concentration, where x = A^dm (7*200 and M Q = 10 13 /i80 _1 M Q is on the order of 
the present-day characteristic mass. The second term in the denominator of equation (21) 
is introduced to prevent the concentration from exceeding a maximum value C max , that 
we set equal to 100, and driving the stellar mass-to-light ratio to very small values. Note 
that B = [C max (l + 7)/C ]~( 1+7 )/ 7 becomes insignificant for 7 <C 1 (e.g, in standard CDM 
models). 

Utilizing the "characteristic density" program kindly provided by J. Navarro, a 
standard cold dark matter (CDM) fluctuation spectrum results in C Q ~ 10 and 7 ~ 0.1. 
For these parameters, r dm /r br ps 150(a viri ai/10) 13 / 30 A~ 5 / 24 (r dm /r e ps 4.5 for a vir i a i rs 10 
and A ~ 1) and r dm is indeed much greater than 10rb r in the luminosity range of interest. 
The dot-dashed curve in Figure 7a shows [M dm /M*] rmax as a function of A for a v i r ; a i = 16 
(i.e. baryon fraction 0.06) CDM models (and for two alternative models to be discussed 
in the next section). Clearly [M dm /M*] rmax increases with optical luminosity, implying a 
decrease of /3 spec with increasing Ly as shown by the dot-dashed curve in Figure 7b. A 
comparison with the observed /3 spec -A correlation (solid curve in Figure 7b), demonstrates 
that a constant baryon fraction CDM model for elliptical galaxy dark matter halos clearly 
fails to match the observed correlation, even over the limited range of A where the 
correlation is well-established. A successful alternative model must produce dark halos 
where [Mdm/M*] rjJUIX decreases with increasing optical luminosity in this range. We will 
discuss the form such possible alternatives must take in the next section. 



5. Discussion 



A natural explanation for why /3 spec , the ratio of the central stellar velocity dispersion 
to the globally-averaged gas temperature, increases with (a) and therefore L v is that the 
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dark-to-luminous mass ratio within the luminous parts of elliptical galaxies (r < r max = 6r e ) 
decreases with L v (see Figures 5, 7a and 7b). In fact, the observed trend requires that 
the total mass-to-light ratio within r max be nearly constant. This can be understood 
from a simple argument based on the virial theorem and FP, as follows. The virial 
temperature within r max , T viria i oc (M dm + M*)/r max oc (l + M dm /M*)(M*/r e ) since r max oc r e 
(all masses are evaluated at r max ). Since, from the FP (a) 2 oc M*/r e , it follows that 
(o-) 2 /T viria i oc (1 + M dm /M*) -1 , or (3 spec oc (1 + M^jM^y 1 if we associate T viria i with the 
integrated X-ray temperature, which should be a fair approximation. 

The implication from the observed I3 spec -Ly trend that the total mass-to-light ratio 
within r max be nearly constant (see below) is contrary to expectations from CDM models 
of large scale structure if all elliptical galaxies have the same global baryon fraction. The 
dark-to-luminous mass ratio is predicted to increase with Ly, in such models, since the 
ratio of the dark matter scale length r dm to the optical half-light radius is large and weakly 
varying, thus producing a predicted correlation between /3 spec and Ly in the opposite sense 
to what is observed (Figures 7a and 7b). We consider two simple variations on the above 
model to try and recover the observed j3 spec -L v trend. First we assume that the dark 
matter concentration depends only weakly on halo mass, as CDM predicts, but relax the 
constant baryon fraction assumption. Then we assume that the concentration varies much 
more steeply with virial mass than the CDM prediction. 



5.1. Alternatives to Constant Baryon Fraction CDM 

If the weak dependence of concentration on halo mass predicted in CDM simulations 
(equation 21 with C a = 10 and 7 = 0.1) is maintained, then the assumption of a constant 
baryon fraction within r2oo must be relaxed: the required increase in dark matter fraction 
within r max for less luminous galaxies must be a reflection of a larger total (i.e. within 
r 2 oo) dark-to-luminous mass ratio. If all galaxies start off with identical baryon fractions, 
this then implies that smaller galaxies lose an increasingly large fraction of their original 
baryonic mass. (Similar conclusions for spiral galaxies have been reached by Persic, Salucci 



|fc Stel 1996 .) We parameterize this variation in the dark/stellar mass ratio as 

(22) 



C^virial C^min 







1 + 







where, as previously defined, atrial = [Md m /M*] r200 , A is the visual luminosity in units of 
L = 5.2 x lO lo /i 8O ~ 2 Ly , and A D is a fitting parameter. Equation (22) has the properties 
^virial — ► 00 as A — > (complete baryonic mass loss), and a V mai ~^ ow n as A — > 00 
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(presumably, (1 + a m i n ) _1 is the average baryon fraction in the universe). The observed 
/3 spec -Ly correlation is satisfactorily reproduced in models with v = 5/4, A Q = 3, and 
«min = 0.5-2.5. The dotted curves in Figures 7a and 7b represent HST stellar models with 
awn = 1-5; the observed trend is very well matched (Figures 7b). (Hernquist stellar models 
require even larger values of v\ however, this is an artifact of the underestimate of f3 spe c for 
low M dm ; see section 4.1). 

Alternatively, if the dark matter halo concentration increases much more steeply 
with virial mass than predicted for a CDM fluctuation spectrum, the observed /3 spec -Ly 
correlation can be reproduced assuming a constant total baryon fraction (i.e., a viria i). In 
this case, the necessary decrease in [M dm /M*] rmax for more luminous galaxies results not 
from a smaller relative global amount of dark matter, but from a larger fraction of the 
dark matter lying outside of the luminous portion of the galaxy (i.e. r max ). We have found 
that equation (21) with 7 = 1, C Q — 4, and a viria i > 5 provides models that adequately 
reproduce the data (smaller values of 7 require larger values of Oviriai)- [M<i m /M*] riiLim and 
{3 S p ec are shown as functions of Ly, for a^mai — 16, by the dashed curves in Figure 7a and 
7b. Such a steep dependence of the concentration on total mass would seem to indicate a 
flatter-than-CDM primordial fluctuation spectrum ( |Navarro, Frenk, fc White 1997| ). 



These alternative scenarios were designed to reproduce the f3 spec -Ly correlation 
without consideration of the scatter or that of galaxies about the mean FP. We address 
these concerns by computing (3 spec for individual pairs (L v , a) corresponding to each of the 
25 non-cluster-cooling-flow, non-point-source-dominated galaxies in the DW sample, for 
the three scenarios of Figures 7a-b. A comparison with the observed values of /3 spec , as 
shown in Figures 8a-c, confirms that the constant baryon fraction CDM scenario provides 
the poorest match to the observations, on average, and further illustrates the magnitude of 
variations about the mean correlation. Further investigation into these variations will be 
pursued in a future work that examines the individual galaxies in detail. 

Both of these alternative scenarios have drawbacks. Mass loss induced by supernovae- 
driven galactic winds during the star-forming epoch of elliptical galaxies has been inferred 
from both the color-magnitude diagram and the enrichment of the intracluster medium 
(e.g., |Loewenstein fc Mushotzky 1996| , pibson 1997|) . However the magnitude of this 



effect must be extreme (i.e u, in equation 22, on the order of 1 or greater) if it is to 
explain the increasing dominance of dark matter in less luminous galaxies required to 
explain the observed T-a relation. The mass loss implied by equation (22) with the 
required set of parameters, when integrated over the elliptical galaxy luminosity function 
( Marzke, Huchra, fc Geller 1994|; |Lin et al. 1996| ), corresponds to ~ 80% of the original 



baryonic mass of L > galaxies. Moreover, in such a scenario, all elliptical galaxies 
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with Ly < 3 x 10 11 /i8oLv Q would span a range in total mass of less than a factor of three; 
the large range in mass loss could not, therefore, be due to less luminous galaxies lying 
within shallower potential wells. Instead, one must presume that variations in the relative 
timing of star formation and the merging of pre-galactic fragments is responsible for the 
large range in mass loss, with galaxies with smaller present-day luminosities having lost the 
largest fraction of their original baryonic mass. Light would be a very poor tracer of mass 
for elliptical galaxies, and the luminosity function would be a reflection of the rarity of 
elliptical galaxies sufficiently relaxed at the star formation epoch to prevent copious mass 
loss from galactic winds. 

It might therefore seem more appealing to resort to a concentration-mass relationship 
for dark halos (equation 21) that is steeper than predicted in CDM simulations in such 
a way that optically less luminous galaxies are more dark matter dominated within r max , 
despite having the same total (within r2oo) baryon fraction. The basic requirement is 
that the dark matter scale length ra m increase from ~ r e at the low-luminosity end of the 
observed range, to r dm ~ r max at L v ~ L = 5.2 x 10 10 /i 8C T 2 Lv/ Q , to r dm ~ r 20 o (C ~ 1) 
at the high luminosity range. These criteria are met for only a rather narrow volume in 
(«viriai, C Q , 7) parameter space, and thus this scenario seems to suffer from a fine-tuning 
problem; nonetheless, the range that successfully reproduces the observed T-cr relation is 
not out of line for cosmogenic models with flat primordial fluctuation spectra ( Navarro" 
Frenk, fc White 19971 ). 



5.2. Mass-to-Light Ratios, Baryon Fractions & Velocity Dispersions 

Figures 9a, 9b, and 9c show the luminosity dependences of mass, mass-to-light ratio, 
and baryon fraction at various radii for the two models discussed above that reproduce the 
observed T-a {f3 spcc -L v ) relation: dotted lines show the model with CDM-predicted halo 
concentration as a function of mass and increasing mass loss with decreasing luminosity 
[C = 10 and 7 = 0.1 in equation 21, a m i n = 1.5, v = 5/4, and A D = 3 in equation 22); 
dashed curves show the model with steeper-than-CDM concentration as a function of mass 
and constant baryon fraction (C = 4 and 7 = 1 in equation 21, a m j n = 16 and A Q = 
in equation 22). We refer to these as "mass loss" and "non-CDM" scenarios, respectively. 
The stellar model is assumed to follow equation (3), and the luminosity range in the plots 
roughly corresponds to that of the DW sample (~ 0.6-2L*). The radii of interest are 
r e = 6.3A 3 / 4 /iso _1 kpc, r^ r = 0.03r e , and r max = 6r e ; r 2 oo^80 increases from 270 to 400 
kpc for the mass loss scenario and r 2 oo ~ 370A 5 / 12 /iso~ 1 kpc for the a v i r i a i = 16, non-CDM 
scenario, where A = L v /L with L = 5.2 x lO lo /i 8O ~ 2 Ly w 3L*. 
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There are a number of robust properties that hold for any model (including Hernquist 
models) to successfully reproduce the T-cr relation. Since the models are constrained 
by the observed properties of the stars and by the average temperature of hot gas in 
hydrostatic equilibrium within r max = 6r e , the inferred quantities at r max are nearly 
model-independent: M(r max )/Lv ~ 25/i8oM /Ly Q , /baryon (?"max) ~ 0.35A 1//4 . Note also that 
M* pa 4.8 lO 11 A 5 / 4 ^8o~ 1 M for any model where the central potential is dominated by the 
stellar component. 

These values agree strikingly well with the constraints from statistical studies of 
gravitational lensing ( |Kochanek 1995| ; [Brainerd, Blandford, fc Smail 1996| ; [Griffiths et al. 



1996). The solid lines in Figure 9a shows the mass as a function of luminosity, evaluated 
at r e and r max for the softened isothermal sphere (SIS) model that provides a best-fit to 
the observed weak gravitational shear from [Griffiths et al. 1996] . It is remarkable that 
the mass-to-light scaling derived from the shear that is dominated by ~ L* galaxies at 
intermediate redshifts is in such excellent agreement with our results using a completely 
independent method on local galaxies primarily with L > L*. For the constant baryon 
fraction CDM model, the total mass within r max is a significantly steeper function of 
luminosity than is demonstrated by the curves in Figure 9a. 

On scales both larger and smaller than r max , the two scenarios depicted in Figures 9a-c 
diverge. In the mass loss scenario the baryon fraction within r e is generally greater than 0.7 
and M{r e ) / L v {r e ) pa 12A 1 / 4 /i 8O M /Lv' q , while in the non-CDM scenario the mass-to-light 
ratio is nearly constant, M(r e )/Ly(r e ) pa ll/i8oM Q /Ly Q , so that the baryon fraction within 
r e becomes less than 0.5 below L pa 0.6L*. 

Extrapolated out to the dark matter virial radius, r 2 oo, the total mass scalings differ 
(by assumption). For the mass loss scenario, 

M virial /L y w 9 [AX 1 / 4 + BX- 1 ] h 80 M Q /L Vc:) , (23) 

where acceptable models have A pa 1.5-5 and B pa 4-8 (A = 2.5, B pa 6 for the model 
shown by dotted curves in Figures 7a-b and 9a-c). For the non-CDM (or any other constant 
baryon fraction) scenario, 

M virial /L y Pa 9O(lO/ baryon )' 1 A 1/4 / i8O M /L y0 ; (24) 

the (constant) total baryon fraction /baryon ~ 0.06 for the dashed curves shown in Figures 
7a-b and 9a-c. We found that /baryon < 0.15 in non-CDM models is required to reproduce 
the observed f3 spec -Ly correlation. Non-CDM models with relatively small global /baryon 
have larger baryon fractions within r e for A < 1 (because of larger dark matter scale lengths) 
and agree somewhat better with the lensing constraints at this radius. 
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Integrating over the early-type galaxy luminosity function flMarzke, Huchra, fc Gellei 



1994| ; |Lin et al. 1996| ) yields an average total (i.e. within the virial radius) mass-to-light ratio 



> 65/i 8 oM /Ly Q for the non-CDM (constant baryon fraction) scenario and fi e ii iptica i s > 0.02, 
which can be compared to an average total mass-to-light ratio of ~ 15O/i8oM0/Ly for 
spirals ( [Persic, Salucci, fc Stel 1996|) , and Spirals = 0.04-0.06. In the mass loss scenario, 



the average total mass-to-light ratio is sensitive to the lower luminosity cutoff, L min , for 
the luminosity function and can be as high as ~ lOOO/i 8 oM©/Ly (^ellipticals ~ 0.4) for 
-^min = 10 8 /i8o~ 2 Ly Q since even the lowest luminosity elliptical galaxy has a total mass of 
~ 3 10 12 /igo _1 M Q . The average density in the universe of virialized mass on galactic scales 
is > 0.07 of the critical density. 

We have calculated the velocity dispersion distributions for both the dark matter and 
stellar density distributions, assuming isotropic orbits. These are compared in Figure 10 
for an Ly = L galaxy for both the mass loss and non-CDM scenarios. Both velocity 
dispersion profiles have maxima since the total gravitational potential is not isothermal. 
The ratio (dark-matter-to-stars) of the squares of these maxima is greater than 1.4 over the 
luminosity range in Figures 9a-c for the two models under consideration, and is ~ 2 over the 
range < Ly < 5L*. In fact, the minimum value of this ratio for any HST or Hernquist 
model that produces /5 spec < 0.7 is greater than one. In this sense, the dark matter is hotter 
than the stars, as predicted by Davis & White (1996). 



6. Summary and Conclusions 

In this work, we have addressed two essential features of the X-ray temperatures 
derived by Davis & White (1996) for an optically complete sample of elliptical galaxies: 
(1) the X-ray emitting gas is always hotter than the stars and, typically, twice as hot 
(Apec = Awip(o") 2 /k(T) ps 0.5); (2) the gas/stellar temperature ratio tends to be higher 
for galaxies with lower velocity dispersions (the "T— a relation"). We have constructed 
physically plausible models of the mass distribution in bright elliptical galaxies in an effort 
to constrain their average dark matter properties by matching these observations. The 
stellar models (described in §3) are fully consistent with the fundamental plane scaling 
relations, and are designed to either conform to the latest published HST results on the 
structure of the centers of elliptical galaxies or to follow the Hernquist (1990) approximation 
to a de Vaucouleurs profile (see Figure 2). We have made other plausible, but conservative, 
assumptions, such as (1) maximizing the non-gravitational (i.e pressure) contribution to 
the gas temperature by allowing the gas and dark matter distributions to extend to infinity 
and (2) assuming that the stellar orbits are isotropic for r <r e . 
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A basic and general result of our calculations is that, in the absence of dark matter, 
Pspec > 0.75 (see Figure 3). Since /9 spcc ~ 0.5 is observed, a convincing case is made that 
dark matter is an extremely common, if not ubiquitous, constituent of elliptical galaxies. 
Although X-ray (and other) observations have been used to infer the presence of dark 
matter in individual cases, we have shown that dark halos are generic to luminous, nearby 
elliptical galaxies. Furthermore, the dark matter /stellar temperature ratios (derived at 
the peak of their velocity dispersion distributions, assuming isotropic orbits) are greater 
than one for models with (3 spec < 0.7. Thus, the observation that the temperature of the 
extended hot gas exceeds the central stellar temperature is a reflection of the fact that the 
dark matter is dynamically "hotter" than the stars, as suggested by Davis & White (1996). 

In §4 we described how /5 spec varies as functions of the two relative parameters of the 
"universal" dark matter distribution (equation 8) — the ratio of dark matter to stellar 
scale lengths, 5 = rdm/?"*, and the ratio of dark-to-luminous matter within r max = 6r e , 
a = [Md m /M*] rmax (see Figures 4a-b and 5). Because we attempt to match only the single 
global observable /3 S pec> there is an allowed range in the details of the dark matter spatial 
distribution; however, we have derived absolute limits on the dark matter parameters 
required to obtain any particular value of /3 spec (see Figures 6a-b). 

The observations do not require that dark matter dominate the inner luminous regions 
of elliptical galaxies: more than half of the mass within r e is baryonic for models with 
Pspec = 0.5 if r dm > r e . 

The most natural explanation of the tendency for galaxies with lower stellar 
temperatures to have larger gas-to-stellar temperature ratios is that [Md m /M*] rmax decreases 
with Ly in such a way that, on average, the total mass-to-light ratio inside r max is 
nearly independent of optical luminosity. This ratio, ~ 25/igoM /Ly , is exactly what is 
predicted for mass models of elliptical galaxies designed to explain the gravitational shear 
of background field galaxies. 

If one specifies a scaling relation for the dark halo concentration, one can extend the 
dark matter distribution out to the virial radius and calculate the total baryon fraction 
and mass-to-light ratio. When we attempt to embed our models within the CDM theory of 
hierarchical halo formation, the implied dark matter scaling badly fails to reproduce the 
observed T-a relation unless smaller galaxies lose an increasingly larger fraction of their 
initial baryonic content (see Figure 7b), such that the average L > galaxy has lost most 
of its initial baryonic mass. Alternatively, the global dark-to-luminous mass ratio could be 
constant if the dark halo concentration declines much more steeply with virial mass than 
CDM models predict, so that the decrease in [Md m /M*] rmax for larger systems is a result of a 
more diffuse dark matter halo rather than a less massive (relative to the stars) halo. In this 
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latter scenario, dark matter may become increasingly important inside r e as L v decreases, 
becoming dominant for L C L». This deviation from CDM predictions of dark halo scaling 
could conceivably be due to a relatively flat primordial fluctuation spectrum on mass 
scales < 10 14 M Q or to the effects on the dark matter density profile of the evolution of the 
baryonic component; to date, large scale structure numerical simulations that can resolve 
halos on galactic scales and include a dissipational component have not been attempted. 

In this paper we have shown that the observed relationship between optical velocity 
dispersions and X-ray temperatures in giant elliptical galaxies implies that they have dark 
matter halos with M/Ly m 25/i8oM /Ly Q within 6r e . In the future, we plan to fully 
utilize the available X-ray and optical imaging and spectroscopic data for our sample on 
a case-by-case basis in order derive the mass distributions in individual galaxies in the 
highest possible detail and to investigate the scatter in the T-a relationship. For galaxies 
with X-ray temperature profiles, we will be able to constrain the detailed form of the dark 
matter distribution as well as its integrated properties. 

We thank Lars Hernquist for reminding us of the consistency considerations of Ciotti 
& Pelligrini (1992), and Richard Mushotzky for feedback on the original manuscript. 
Comments from an anonymous referee led to significant improvements in the quality of this 
paper. R.E.W. acknowledges partial support from NASA grants NAG 5-1718 and NAG 
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A. Projection of a Broken Power-law 



Suppose the stellar density profile is given by a broken power-law of the form 




(Al) 



The projected surface mass, 




(A2) 



can be recast, utilizing the change of variable y = (r 2 — R 2 )/R 2 , into the form 



£*(&) 



p* G r c J Ib-^-Vil+b)-* , 



(A3) 



where 




(A4) 
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b = R 2 /r 2 , and c = 6(1 + b)^ 1 . Note that the bracketed term in equation (A3) has the 
asymptotic slopes at small and large R that one would expect from dimensional analysis 
(i.e. £* oc p*r). The definite integral of equation (A4) can be evaluated as follows 
( Pradshteyn and Kyzhik 1980 ): 



J = c ~" B( 7 _ 0.5, 0.5) 2 F 1 (a, 7 - 0.5; 7; -6" 1 ), (A5) 
where 7 = 0" + 0.5k, is one-half the asymptotic density slope as r — > 00, 5 is the beta 



function, and 2 -P\ is the hypergeometric function (|Gradshteyn and Ryzhik 1980Q . We have 



empirically found that we can introduce a parameter r, recast equation (A3) as 

= p* r c J'b-^- T \l + b )-^-°- 5 \ (A6) 

where 

and determine a value of r (for any k < 1) that results in a very nearly constant J' as a 
function of R (with no significant structure near R = r c ). Equation (A6) has the expected 
slope as R — > 00; and, we find r — > as k — >• (the "/3 model"), and r — > 0.5 as k — >• 1. We 
conclude that the projection of a density profile of the form expressed by equation (Al) 
is well-described by equation (2) (at least for a = 2), thus perhaps lending some physical 
justification to its use. 
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Fig. 1. — (3 spec vs. log L v for the DW sample (points, with vertical lines representing 
90% confidence limits), and best-fit /5 spec - log L v (dashed curve) and log(/5 spec ) - log L v 
(dot-dashed curve) linear regressions. See Section 4.5 for details. 

Fig. 2. — Comparison of luminosity density profiles, Z(r), for "deV" (solid curve), "HST" 
(dotted curve), and "Hernquist" (dashed curve) stellar models (see text). The radii are 
normalized to r max = 6r e , l(r) to the average luminosity density within r max , l avg = 
3L/47rr max 3 where L is the total stellar luminosity. 

Fig. 3. — p s p ec vs velocity dispersion anisotropy scale length, s, in units of the break radius, 
rbr, for HST (a = 0: solid curve, a = 3: dashed curve), and Hernquist (a = 0: dotted curve, 
a = 3: dot-dashed curve) stellar models, (a = [M dm /M^] rmax .) 

Fig. 4. — (a) (3 spec vs r^ra/r^ for [Md m = 3M*] rmax (solid curve: HST stellar model, dotted 
curve: Hernquist model), (b) Baryon fractions at rt, r (dashed curve for HST, dot-dashed 
curve for Hernquist model) and r e (solid curve for HST, dotted curve for Hernquist model) 

vs r dm /r br for [M dm = 3M*] rmax . 

Fig. 5. — p s p ec vs [Md m /M*] rmax for HST (solid curves), and Hernquist (dotted curves) stellar 
models and r d m =0.5, 1.0, and 2.0r e (upper, middle, and lower curves, respectively). 

Fig. 6. — (a) Minimum values of rdm/^br for HST and Hernquist models, (b) Maximum values 
of baryon fractions within (from top to bottom) r br , r e , and r max for HST and Hernquist 
models. The HST (Hernquist) model is represented by solid (dotted) curves. 

Fig. 7. — (a) Ratio of dark to luminous matter within r max , vs dimensionless luminosity, 
A = L v /L (L = 5.2xl0 w h 80 - 2 L Vl3 ), for HST stellar models. The dot-dashed curve denotes 
constant (within the virial radius) baryon fraction (/baryon = 0.06) model and CDM scaling 
of dark matter concentration; dotted curve has /baryon increasing with optical luminosity; 
dashed curve has a steeper-than-CDM scaling of concentration with dark halo mass, (b) 
Predicted variation of (3 spcc with dimensionless luminosity. Line-types are as in (a). Data 
points and best-fit correlation (solid curve) from Figure 1 are re-plotted. (For clarity, three 
galaxies where f3 spec is uncertain to > 50% are omitted from this and the following figure.) 

Fig. 8. — (a) /3 spcc calculated using pairs (Ly, a) as measured for the DW sample and CDM 
scaling of dark matter concentration with constant /baryon vs. observed [3 spec . (b) Same as 
(a) for model where /baryon increases with optical luminosity, (c) Same as (a) for model with 
steeper-than-CDM scaling of concentration with dark halo mass. The solid line represents 

spec fiobs- 



-27- 



Fig. 9. — (a) Mass vs L v (in solar units) within, from bottom to top, r e , r max = 6r e (stars 
only), r max (total), and r 2 oo- The solid curves show the mass at r e and r max from weak 
gravitational lensing studies, (b) Same as (a) for mass-to-light ratio (in solar units) at (from 
lowest to uppermost curves): r = 0, r e , r max , and r2oo- ( c ) Same as (a) for baryon fraction 
vs Ly at (from uppermost to lowest curves) r = r^ T , r e , r max , and r 2 oo- Dotted and dashed 
line-types have the same connotations as in Fig. 7. 

Fig. 10. — 1-d velocity dispersion distributions, assuming isotropic orbits, for an L v = L = 
5.2 x 10 10 /i 8 o _2 Ly Q galaxy. "Mass loss" and "non-CDM" models are denoted by dotted 
and solid curves, respectively, for the stellar profiles and dot-dashed and dashed curves, 
respectively, for the dark matter profiles. [Md m /M*] rmax = 1.5(1.3) and r<i m /r e = 4.1(12) 
for the mass loss (non-CDM) models at this luminosity, and the (observational and model) 
value of P spcc is 0.62. 
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